Steady-state Levy flights in a confined domain 
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Abstract 

We derive the generalized Fokker-Planck equation associated with a Langevin equation driven 
by arbitrary additive white noise. We apply our result to study the distribution of symmetric and 
asymmetric Levy flights in an infinitely deep potential well. The fractional Fokker-Planck equation 
for Levy flights is derived and solved analytically in the steady state. It is shown that Levy flights 
are distributed according to the beta distribution, whose probability density becomes singular at 
the boundaries of the well. The origin of the preferred concentration of flying objects near the 
boundaries in nonequilibrium systems is clarified. 
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I. INTRODUCTION 



It is remarkable, and at first sight surprising, that a large variety of physical, biological, 
: inancial, and other processes can be described by stable Levy processes with infinite variance 
The latter are defined as continuous-time random processes whose independent 
and stationary increments are distributed according to heavy-tailed stable distributions. 
The main feature of these distributions is that the tails cannot be cut off or, in other words, 
rare but large events cannot be neglected. As a consequence, the classical stochastic theory, 
which is based on the ordinary central limit theorem, is no longer valid. 

Due to the heavy-tailed distributions of the increments, stable Levy processes exhibit 
large jumps, and for this reason these processes are often called Levy flights. Levy flights 
are actually observed in various real systems. Representative examples include, for instance, 



fluorescent probes in living polymers 4j, tracer particles in rotating 
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lattices jg], cooled atoms in laser fields 7], subsurface hydrology 



lows [5J , ions in optical 



9] 



and ecology 



id, 



12j , though a recent work questions the empirical evidence for Levy flights in animal 
search patterns Levy flights have also been predicted for a large number of model 

systems The ubiquity of these processes is supported by the generalized central 

limit theorem |l4j, which states that all limiting distributions of properly normalized and 
centered sums of independent, identically distributed random variables are stable. 

The Langevin equation is one of the most important tools for studying noise phenomena 
in systems coupled to a fluctuating environment. Introduced by Paul Langevin just one 



hundred years ago 15) to describe the dynamics of a Brownian particle, this equation and its 
various modifications are widely used in many areas of science [l6]. The Langevin equation 
driven by Levy white noise, i.e., noise defined as a time derivative, in the sense of generalized 
functions, of a stable Levy process, provides a basis for the study of Levy flights in external 



potentials. It has been shown that t he probab i 
fractional Fokker-Planck (FP) equation 
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ity density of Levy flights satisfies the 



20, 



2l|. The steady-state solutions of this 



equation describing confined Levy flights, i.e., flights with finite variance, are of particular 
interest. One reason is that these solutions will clarify the distribution of flying objects in 
confined domains. This is an important issue, especially near impermeable boundaries, in 
such complex systems as confined plasmas and turbulent flows. Another reason is that exact 
general solutions of a simple form, which are valid for any Levy white noise, will be very 
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useful for testing a variety of numerical methods in this area 22|, |23|, |24|. However, known 



solutions are related to power potentials and to a very special case of Levy white noise with 



unit index of stability and zero skewness parameter [25|, |26| , and thus they are not suitable 
for those purposes. 

It should be noted that Levy white noises do not exhaust all possible white noises. As 
a consequence, the fractional FP equation is a particular case of the generalized FP equa- 
tion, which corresponds to the Langevin equation driven by an arbitrary white noise. Since 
any white noise is defined as a time derivative, in the sense of generalized functions, of a 
stationary process with independent increments, it can be characterized by the transition 
probability density or, alternatively, by the characteristic function of this white noise gen- 
erating process. One expects therefore that the term in the generalized FP equation that 
describes the effect of the noise on the dynamics of the system can also be expressed via the 
characteristic function. The derivation of the generalized FP equation is of great importance 
because it accounts for all possible white noise effects in a unified way and will be very useful 
for applications. 

In this paper, we put forward the generalized FP equation and find the analytical solution 
of this equation for steady-state Levy flights in a confined geometry. 

II. GENERALIZED FOKKER-PLANCK EQUATION 

In many applications, ranging from physical and chemical to biological and social systems, 
the relevant degrees of freedom of these systems obey a (dimensionless) Langevin equation 
that is equivalent to the equation of motion for an overdamped particle 

x(t) = f(x{t),t)+£(t). (1) 

Here, x(t) is a particle coordinate, x(0) = 0, f(x,t) = —dU(x,t)/dx is a force field, U(x,t) 
is an external deterministic potential, and £(£) is a random force (noise) resulting from 
a fluctuating environment. Though the quantities in Eq. (TjQ) have different meanings for 
different systems, we will use the above terminology to be concrete. 

Under certain conditions (see, e.g., Refs. 0, |28j), the noise £(t) can be chosen to be 



white. In this case, the increment Sx(t) = x(t + r) — x(t) of the particle coordinate during 
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a time interval r (r — > 0) is written as 

fa(t)=/(s(t),t)T + <fy(t), (2) 

which defines the meaning of Eq. ([1]) in the white- noise approximation. Here Srj(t) = 
rj{t + r) — rj{t) = f^ +T dt'£(t'), and we assume that the integral exists in the mean 
square sense. A white noise generating process, i.e., a stationary random process rj(t) = 
lim T ^o Y^}j=o~ l ^(j r )) with 77(0) = and where [t/r] denotes the integer part of t/r, is 
completely defined by the transition probability density p(rjj + i, T\rjj) of a discrete-time pro- 
cess rj{nr) = Y^jZo dv(j T ) (n > 1) as t ^ 0. Here r)j + i and r]j denote the possible values of 
r)(jr + r) and t)(Jt), respectively. Note that all transition probability densities of the form 



equation 



77,-) can be expressed through p(r)j + i,T\r)j) by using the Chapman- Kolmogorov 
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3l|. In particular, p^^^r^) = J dr) j+ ip(r) j+2 , T\r) j+1 )p(r] j+ i, r|^). 
This implies that the influence of any white noise on the system can also be fully char- 
acterized by the function p(rij + i,r\rjj). We assume that the transition probability density 
p(rjj +1 ,T\rjj) is properly normalized, J™ dr)j + ip(r]j + i, r\r]j) = 1, and that it satisfies the con- 
dition lim T ^Q p(rjj + i, r\rjj) = — r)j), where 5(-) stands for the Dirac 5 function. For 
simplicity, we also assume that p(r)j + i, r\r]j) = p(Arj, r) with Ar] = r]j + i — rjj. 

We define the probability density of the particle coordinate x(t) in the usual way, 
namely P(x,t) = (5(x — x(t))), where the angular brackets denote averaging over the 
noise. Taking the Fourier transform of P(x,t) according to the definition T{g{x)} = g k = 
dx e~ tkx g(x), we obtain P k (t) = (e~ lfca:( -^), i.e., the characteristic function of x(t). Equa- 
tion d2J) implies that the increment of this quantity, 5Pk = P k (t + r) — P k (t), can be writ- 
ten in the form 5P k = -ikr{e- ikx{ $ f{x{t), t)) + (e~ ikx ^ ( e ~ ikS ^ - 1)) as r -> 0. The 
use of the well-known properties of the Fourier transform yields ik(e~ tkx ^ f(x(t),t)) 



J 7 {df(x,t)P(x,t)/dx}, and the statistical independence of x(t) and Sr](t) implies [29 ] 
( e -ikx(t) ( e -iks v (t) _ ^ = p k (t)(j) k (r) - 1), where p k (r) = T{p{At),t)} = ( e - ik5r >^) is the 
characteristic function of Sr](t). Dividing 5P k by r and taking the limit r — > 0, we obtain 
the generalized FP equation in Fourier space, 

-P k {t) = -F[—f(x,t)P(x,t)} + P k {t)ct> k , (3) 

where fc = lim^oO^O) - l)/r. 

It is advantageous to introduce the characteristic function S k = (e~ tkv ^) of r/(l). We 
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rewrite it as S*. = lim^oQo^r))' 1 ^, using the formula 7/(1) = lim r ^ X^=o' _1 ^v(j T )- R- e_ 
placing pjfc(r) by 1 + rcpk and taking into account that lim e ^ (l + £) 1//£ = e, we readily 
find that Sk = e™. Finally, applying the inverse Fourier transform T^ x {gk\ = g{x) = 
(l/2?r) JZodke^gk to Eq. © and using 0^ = lnSfc, we obtain the generalized FP equation 
in real space 

(4) 



-P(M) = -— f(x,t)P(x,t) +F- 1 {P k (t)lnS k }, 



with P(x, 0) = <5(x), which corresponds to the Langevin equation (0Q) driven by an arbitrary 
white noise. 

Equation (j3J) represents our first main result. It constitutes a closed, concise representa- 
tion of the combination of the Fokker-Planck and Kolmogorov-Feller equations, which are 



the 



21, 



3asic equations governing continuous and discontinuous Markov processes, respectively 
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3l|. A remarkable feature of this equation is that it accounts for the noise influence 



in a unified way, namely by means of the characteristic function Sk of the white noise gen- 
erating process rj(t) at t = 1. All presently known FP equations associated with Eq. (pQ) can 
be obtained directly from Eq. (j3J). In particular, if £(£) is Poisson white noise characterized 
by the transition probability density p(Arj, r) = (1 — Ar) S(Arj) + Arg(A?7), where A is the 



average number of jumps of rj(t) per unit time and 



sizes, then Sk = e A ^ 9fe ) and Eq. 





Of 



P(x,t) 



yields 
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[Ar}) is the probability density of jump 



34] 



dx 



+A 



f{x,t)P{x,t) - XP(x,t) 



dyP(y,t)q(x-y). 



(5) 



If £ (t) is Levy white noise then the generalized central limit theorem |l4j implies that Sh is 
the characteristic function of Levy stable distributions. As is well known (see, e.g., Ref. [35]), 
the characteristic function Sk{a, j3, 7, p) of non-degenerate stable distributions depends on 
four parameters: an index of stability a G (0,2], a skewness parameter (3 G [—1, 1], a scale 
parameter 7 G (0, 00), and a location parameter p G (—00,00). Assuming in accordance 
with the initial condition P(x, 0) = 5(x) that p = and excluding from consideration the 



singular case where a = 1 and (3 7^ simultaneously (in this case \(f>k 
reaches the final state immediately), we obtain Sk = Sk («,/?, 7), where [351] 



00, and the system 



S k (a,(3,j) = exp 



-7|fc| a ( 1 + i(3 sgn(k) tan — 



(6) 
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Equation (j4]) with = Sfc(a,/3,7) can be easily rewritten as a fractional differential 
equation. The Riemann-Liouville derivatives of a function g(x) on the interval [— s, s] are 



defined as 



36] 



,D±g(x) 



f±l) n d n 



dyg(xTy) V 



n—cr—l 



(7) 



T(n — a) dx jo 

Here, S D^_ and S D°_ denote the operators of the left- and right-hand side derivatives of order 
o (0 < o < oo), respectively, with n — 1 + [a], and T(z) is the gamma function. Since 
^{ooD^P^x, t)} = (±ik) a Pk(t), which follows from the definition ([70, Eq. (jlj) reduces to the 
desired fractional FP equation: 



d_ 

dt 



P(x,t) 



~f(x,t)P(x,t)-- j- 

ox 2cos(7tq;/2) 



'l-P) 00 D1}P(x,t). 



(8) 



All previously known forms of the fractional FP equation, which correspond to the Langevin 
equation ([T|) driven by Levy white noise, can be derived from Eq. ([8]). In particular, taking 
into account the relation (oo-D? + ^D^Pix, t) = 2 cos(7ra/2) J 7 ' 1 {\k\ a P k {t)} and the defi- 



nition of the fractional Riesz derivative 



36j, d a P(x,t)/d\x 



the case of symmetric Levy white noise {(3 = 0) yields 17 
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2l| 



a P k (t)}, Eq. © in 



§£P(X, t) = -±f( X , t)P(x, t) + 1^P(X, t). 



(9) 



Specifically, if represents Gaussian white noise of intensity D, i.e., a = 2 and 7 = D, 
then Eq. <Q becomes the ordinary FP equation (37) ] . 

We note that Eq. (jlj) also applies to the case of compound white noises. For example, if 
£(t) is the sum of independent Levy and Poisson white noises then Sk = Sk(a, /3, 7) e~ x ^~ qk ^ 
and Eq. (jlj) can be written as 



0_ 

Of 



P(x, t) = ~f{x, t)P(x, t) - ^— - [(1 + (3) 

ox 2cos(TTa/ 2) 

x 00 D^ + (l-f3) 00 D a _}P(x,t) 

POD 

-XP(x,t)+X dyP(y,t)q(x-y). 



(10) 



III. STEADY-STATE LEVY FLIGHTS IN A CONFINED GEOMETRY 



We apply Eq. (jSJ) to the case of stationary Levy flights in an infinitely deep potential well. 
We assume that f(x,t) = within the well, i.e., for x G [—1,1], and that the boundaries at 
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x = ±Z are impermeable for particles, i.e., P(x,t) = at \x\ > I. With these conditions, 
Eq. (jSJ) for the stationary probability density P s t(x) reduces to (1 + j3)iD"P st (x) + (1 — 
j3)iD"P st (x) = 0. Rewriting this equation as dJ(x)/dx = 0, where J(x) is the probability 
current, and using the boundary condition J(±Z) = [37], we obtain the equation J(x) = 0, 
which for < a < 1 reads 



The fact that iD±(l ± = [36( suggests seeking a solution of Eq. (TTT1) in the form 

-fst(^) — C(Z + x)~ u (l — where C is a normalization factor. The parameters v and /x 

are determined by the equation 

{l+(3)B(l-a,l-u)z l - a F(l-a,fi;2-a-u; -z) 

= (l-/3)B(l-a, l-^)F(l-a, i/; 2-a-/i; (12) 

Here, z = (I + x)/(l — x), B(a, 6) = r(a)T(6)/r(a + 6) is the beta function, and F(a, b; c; y) 
is the Gauss hypergeometric function. Equation ffl2l must be independent of x, since v and 
fi do not depend on x. This requirement leads to the condition a + v + /x = 2. Using the 
relation F(a, 6; 6; — y) = (1 + y)~ a [38], we find that in this case Eq. ( fl2"l) becomes 

(1 + /3)B(1 - a, 1 - 1/) - (1 - /3)B(1 - a, 1 - /i) = 0. (13) 

Solving the equations a + z/ + ji = 2 and (TT3"]) with respect to /i and 1/, we find 

I I = 1 — — ± — arctan \J3 tan — J , (14) 

where arctan x denotes the principal value of the inverse tangent function. Similar calcu- 
lations for 1 < a < 2 lead to the same result, and formula (|14p is valid for all a £ (0, 2] 
(excluding the case a = 1, (3 ^ 0). Finally, calculating the normalization factor C, we obtain 

P«W = ( 2l ) 1 ~ a il nn ) ~ U{ \~ X) 7 - (15) 

S(l -1/,1-fi) 

Equations (|14p and ffTol) represent our second main result. Particles that perform Levy 
flights are distributed in an infinitely deep well according to the beta distribution (see Fig. 1). 
The main feature of this distribution is the singular behavior of P 8 t(x) as |x| — > I if a < 2 
and (3 7^ ±1. The reason is that for a < 2 the particles can perform random jumps in both 



directions. However, the boundaries are impermeable, and consequently the particles are 
concentrated preferably near these two boundaries. In particular, for (3 = 0, Eqs. (TH1) and 
(USD yield P st (x) = (2l) 1 -° i T(a){l 2 - x 2 )^ 2 " 1 /V 2 (a/2). In contrast, for (5 = ±1, one-sided 
jumps dominate and the particles are concentrated near one of the boundaries. Specifically, 
P st (x) = 5(1 - sgn(p)x), if < a < 1, and P st (x) = (2f) 1_af (a - 1)(7 + sgn(p)x) a ~ 2 , if 
1 < a < 2. Finally, for a = 2 the sample paths of x(t) are continuous and the stationary 
distribution is uniform, i.e., P s t(x) = 1/(2/). 

IV. CONCLUSIONS 

We have derived the generalized FP equation associated with the Langevin equation 
driven by an arbitrary white noise. This FP equation accounts for the influence of the noise 
by means of the characteristic function of the white noise generating process. In the case 
of Levy flights, this equation has been reduced to a fractional FP equation and has been 
solved analytically in the steady state for a confined domain. It has been shown that both 
symmetric and asymmetric Levy flights in an infinitely deep potential well are distributed 
according to the beta probability density. The preferred concentration of flying objects near 
impenetrable boundaries results from the jumping character of Levy flights. 
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FIG. 1: (Color online) Plots of the stationary probability density (|15p for different values of the 
parameters a and f3 and 1 = 1. 
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